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^ ■ Abstract. We investigate spectral properties of a quantum particle confined to an infinite 

O \ straight planar strip by imposing Robin boundary conditions with variable coupling. As- 

■ suming that the coupling function tends to a constant at infinity, we localize the essential 
spectrum and derive a sufficient condition which guarantees the existence of bound states. 

■ Further properties of the associated eigenvalues and eigenfunctions are studied numerically 
by the mode-matching technique. 
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S: 

1 Introduction 

\ Modern experimental techniques make it possible to fabricate tiny semiconductor structures 

I which are small enough to exhibit quantum effects. These systems are sometimes called nano- 

■ structures because of their typical size in a direction and they are expected to become the 
building elements of the next-generation electronics. Since the used materials are very pure 
and of crystallic structure, the particle motion inside a nanostructure can be modeled by a free 

■ particle with an effective mass m* living in a spatial region Vl. That is, the quantum Hamiltonian 
. can be identified with the operator 



2m* 

in the Hilbert space L^(r2), where h denotes the Planck constant. We refer to |1H [23] for more 
information on the physical background. 

An important category of nanostructures is represented by quantum waveguides, which are 
modeled by Q being an infinitely stretched tubular region in or M'^. In principle, one can 
consider various conditions on the boundary of in order to model the fact that the particle 
is confined to 17. However, since the particle wavefunctions ip are observed to be suppressed 
near the interface between two different semiconductor materials, one usually imposes Dirichlet 
boundary conditions, i.e. ip = on dil. Such models were extensively studied. The simplest 
possible system is a straight tube. The spectral properties of corresponding Hamiltonian in this 
case are trivial in the sense that the discrete spectrum is empty. 

It is known, that a deviation from the straight tube can give rise to non-trivial spectral 
properties like existence of bound states, by bending it [U [HI [111 [TBI EH [22l [24] , introducing an 



*This paper is a contribution to the Proceedings of the 3-rd Microconference "Analytic and Algebraic Me- 
thods III" . The full collection is available at |http: / /www. emis.de /journals/SIGM A /Prague2007.html , 
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Figure 1. Straight waveguide with Robin boundary conditions. The dotted hnes indicate that the 
boundary-couphng function is aUowcd to vary. 

arbitrarily small 'bump' [3l [5] or impurities modeled by Dirac interaction [l^, coupling several 
waveguides by a window [16], etc. 

Another possibility of generating bound states is the changing of boundary conditions. It can 
be done by imposing a combination of Dirichlet and Neumann boundary conditions on different 
parts of the boundary. Such models were studied in [9l [101 [13 ED]- 

In this paper, we introduce and study the model of straight planar quantum waveguide with 
Robin conditions on the boundary. While to impose the Dirichlet boundary conditions means 
to require the vanishing of wavefunction on the boundary of 0, the Robin conditions correspond 
to the weaker requirement of vanishing of the probability current, in the sense that its normal 
component vanishes on the boundary, i.e. 

j • n = on di}, 

where the probability current j is defined by 

This less restrictive requirement may in principle model different types of interface in mate- 
rials, namely, a superconducting film sandwiched between two metallic non-superconducting 
substrates. 

The Laplacian subject to Robin boundary conditions can be also used in the problem to 
find the electro-magnetic field outside the object consisting of a conducting core covered by 
a dielectric layer. If the thickness of the layer is too small compared to the dimension of the 
conductive core, the numerical methods of solving this problem fail because of instabilities that 
then arise. In this case, the problem can by solved by approximation of the dielectric layer by 
appropriate boundary conditions of the Robin type. We refer to [21 [12] for more information. 

The system we are going to study is sketched on Fig. [T) We consider the quantum particle 
whose motion is confined to a straight planar strip of width d. We shall denote this configuration 
space by := M X (0, d). For definiteness we assume that it is placed to the upper side of the 
X-axis. On the boundary the Robin conditions are imposed. More precisely, we suppose that 
every wave-function ip satisfies 

—dyip{x, 0) -I- a{x)'il;{x, 0) = 0, 

dyi^ix, d) + a{x)ij{x, d) = 0, (1.2) 

for all x G M. Notice that the parameter a depends on the x-coordinate and this dependence is 
the same on both "sides" of the strip. We require that a{x) is positive for all x S M. Moreover, 
in Theorem [1] we will show that a sufficient condition for the self-adjointness of the Hamiltonian 
is the requirement that a G Vl^^'°°(R). 

Putting h'^/2m* = 1 in ([TH), we may identify the particle Hamiltonian with the self- adjoint 
operator on the Hilbert space L^(J7), defined in the following way 

HaTp := -AtP, VV' G D{Ha) := G W'^''^{n) \ ^ satisfies O} , (1-3) 

where D(Ha) denotes the domain of the Hamiltonian. 
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While it is easy to see that Ha is symmetric, it is quite difficult to prove that it is self-adjoint. 
This will be done in the next section. SectionOis devoted to localization of the essential spectrum 
and proving the existence of the discrete spectrum. In the final section we study an example 
numerically to illustrate the spectral properties. 



2 The self-adjointness of the Robin Laplacian 

For showing the self-adjointness of the Hamiltonian, we were inspired by [H Section 3]. Our 
strategy is to show that is, in fact, equal to another operator H^, which is self-adjoint and 
defined in following way. 

Let us introduce a sesquilinear form 



ha{4i,il^) ■= / V(/)(x,?/) • V^(x,y) dxdy -I- / a{x)(^(f)(x,0)'ip{x,0) + (f){x,d)ip{x,d)) dx 
Jn Jr 

with the domain 

D{ha) := W^^^in). 

Here the dot denotes the scalar product in and the boundary terms should be understood 
in the sense of traces [H Section 4]. We shall denote the corresponding quadratic form by 
•= ha{ip, tp). In view of the first representation theorem |19l Theorem VI. 2.1], there exists 
the unique self-adjoint operator Ha in L'^{^) such that haifj), V') = Ha'4') for all ^ € D{Ha) C 
D{ha) and (f) G D{ha), where 

D{Ha) = G D{ha) I 3F G L\n),ycj) G D{ha),ha{(l),ip) = {(p,F)} . 

For showing the equality between H^ and H^ we will need the following result. 

Lemma 1. Let a G W^'°°{U.) and Vx G M, a{x) > 0. For each F G L'^{^), a solution ip to the 
problem 

ycpGW^'\n), ha{(t>,^) = {<t),F) (2.1) 

belongs to D{Ha)- 

Proof. For any function tp G W^''^{Q,), we introduce the difference quotient 

ip{x + 5,y) -ilj{x,y) 
ips{x,y) ■■= , 

where (5 is a small real number. Since 



\i;{x + 6,y) - ij{x,y)\ 
we get the estimate 



^ / dxip{x + 5t, y) dt 







< \6\ / \dxij{x + 6t,y)\ dt, 







f |^5|'< / { C \dxi^{x + 5t,y)\ dt] dxdy< I ( T \dxi;{x + 5t,y)f dt ] dxdy 
Jn Jn \Jo / Jn \Jo 

= [' ( [ \d.ijix + 5t,y)f dxdy) dt = [ |5,V'(x,y)|2 drcdy. 
Jo \Jn / Jn 
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Therefore the inequahty 

Ush-iu) < UWw^^-in) (2.2) 
holds true. 

If ip satisfies (j2.ip . then is a solution to the problem 



ha{4>, i^s) = {4>, Fs) - / as{x) I (/)(x, 0)'ip{x + 6,0) + 4>{x, d)'4}{x + 6,d) \ dx, 

where <p G VF^'^(r2) is arbitrary. Letting (p = Tps and using the "integration-by-parts" formula 
for the difference quotients, {(j),Fs) = — (<^_5,F), we get 



K[iJ5] = -{{il^s)-S,F) - / as{x) {i;s{x,0)i;{x + 6,0) + M^,d)i;{x + 6,d)) dx. (2.3) 



Using Schwarz inequality, Cauchy inequality, estimate (j2.2p . boundedness of a and as, and 
embedding of VF^'^(r2) in L'^{d^l), we can make following estimates 



1—2 , 1,,.,, „2 



\{{ijs)^5,F)\ < \\F\\L2^n)\\{'ilJ5)-5\\LHn) < + i^Us\\w^,2(^q), 



as{x) yiljs{x, 0)'ip{x + 6,0) + ips{x, d)ip{x + 6, d)j dx 
< Ci\\'>ps\\L2(dn)\\4'\\L2(dn) < C'2||V'5|liyi.2(n)||'0llvyi.2(n), 
{x){\M^,0)\^ + \Mx,d)\^) dx 



a 



< C's 11^5 1 1^2 (9!^) < C4||V'5||^i,2(n) 

with constants C1-C4 independent of 6. Giving this estimates together, the identity (|2.3p yields 

-||V'5||^i,2(Q) < C4||V'5||^i.2(Q) +(^2 11^/^5 |ltyi,2(o)llV'llH^i>2(o) + 2ll-^llL2(f7) + 2^\Mwi.2(n)- 
We get the inequality 

where the constant C is independent of 6. This estimate implies 
sup ||V'-5||H/i.2(m < 00; 

and, therefore, by \13\ § D.4] there exists a function v € VK^'^(r2) and a subsequence (5^ — > 
such that '(/'-^j. ^ I) in W^''^{il.). But then 

- / dxipcp = / ^^xf^ = / V lim = lim / = - lim / V-5</' = - / t^<A- 

Jn JQ JQ Ofc-»0 Ok^OjQ Sk^OjQ JQ 

Thus, dxip = V in the weak sense, and so dxip G W^''^{Q). Hence, dxx'fp £ L^(r2) and f?!,^';/' G 
L2(f7). 

It follows from the standard elliptic regularity theorems (see [131 § 6.3]) that ^ G TVj^'^^(r2). 
Hence, the equation —Aip = F holds true a.e. in Thus, dyyip = —F — dxxip £ F'^i^): and 
therefore ip G VF^'^(r2). 

It remains to check boundary conditions for ^. Using integration by parts, one has 



(0, F) = /i„(^,0) = (,/.,- AV')+ / (t>{x,0)[-dy'^{x,0) +a{x)i;{x,0)] dx 

Jr 

+ / (j){x,d) [dyi/j{x,d) + a{x)'il;{x,d)] dx 
Jr 

for any cp G VF^'^(r2). This implies the boundary conditions because —Aip = F a.e. in Q and 
is arbitrary. I 
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Theorem 1. Let a G W^^°°{'K) and^x e M, a{x) > 0. Then = Ha- 

Proof. Let -0 € D{Ha), i.e., G W'^''^{fl) and ^ satisfies the boundary conditions (|1.2|) . Then 
ip G D{ha) = W^''^{Q) and by integration by parts and (|1.2p we get for all cp G D{ha) the 
relation 



ha{4>,ip) = / (t>{x,d)dyil){x,d) Ax — / (j){x,Q)dy'4>{x,Q) dx — / (^( 

Jn 

+ I a{x)(f>{x,d)'ip{x,d) dx + / a{x)(l){x,0)'ilj{x,0) dx 



x,y)Aip{x,y) dx dy 



(j){x, y)Ail;{x, y) dx dy. 

It means that there exists rj := —Ai/j G such that V0 G D{ha), ha{4>,ip) = {(j^^v)- That 

is, is an extension of Ho,. 

The other inclusion holds as a direct consequence of Lemma [T] and the first representation 
theorem. ■ 



3 The spectrum of Hamiltonian 

In this section we will investigate the spectrum of the Hamiltonian with respect to the behavior 
of the function a. In whole section we suppose that a G W^'°^{R) and Vx G M, a{x) > 0. We 
start with the simplest case. 



3.1 Unperturbed system 

If a{x) = ao > is a constant function, the Schrodinger equation can be easily solved by 
separation of variables. The spectrum of the Hamiltonian is then a{HaQ) = [-Bi(ao)) oo), where 
i?i(ao) is the first transversal eigenvalue. The transversal eigenfunctions have the form 

a 



Xn{y;a) = Na y—^== sin (^^/Ej^yj + cos (^y/ Enia)yj j , (3.1) 

where is a normalization constant and the eigenvalues En{a) are determined by the implicit 
equation 



f{En;a) = 2a^/En{a) cos{^/En{a)d) + (a^ - ^„(a)) sin( V^„(a)d) = 0. 
Note that there are no eigenvalues below the bottom of the essential spectrum, i.e., ddisc = 0- 



3.2 The stability of essential spectrum 

As we have seen, if a is a constant function, the essential spectrum of the Hamiltonian is the 
interval [-Ei(ao), oo). Now we prove that the same spectral result holds if a tends to oq at 
infinity. 

Theorem 2. // lim a{x) — oq = then (Tcss{Ha) = [Ei(ao), oo). 

\x\^oo 

The proof of this theorem is achieved in two steps. Firstly, in LemmaO we employ a Neumann 
bracketing argument to show that the threshold of essential spectrum does not descent below 
the energy Ei{ao). Secondly, in LemmaHJ we prove that all values above Ei{ao) belongs to the 
essential spectrum by means of the following characterization of the essential spectrum which 
we have adopted from [8|. 
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Lemma 2. Let H be a non-negative self-adjoint operator in a complex Hilbert space 7i and h 
he the associated quadratic form. Then A € acss{H) if and only if there exists a sequence 
{'4'n}'^=i C D{h) such that 

(i) V?iGN\{0}, llVnII = 1, 
(a) ilJn — - — > in TL, 

n—*oo 

(Hi) {H - \)ipn > in D{hY . 

n— >oo 

Here D[h)* denotes the dual of the space D{h). We note that 

II /II K'^.V')! 
WmiDih)* = sup II II 
<t,£D{h)\{o} \m\i 

with 

The main advantage of Lemma [2] is that it requires to find a sequence from the form domain 
of H only, and not from D{H) as it is required by the Weyl criterion [71 Lemma 4.1.2]. Moreover, 
in order to check the limit from (Hi), it is still sufficient to consider the operator H in the form 
sense, i.e. we will not need to assume that a is differentiable in our case. 

Lemma 3. If lim a{x) — qq = then ini acss{Ha) ^ Ei{cto)- 

\x\~*OD 

Proof. Since a{x) — ao vanishes at infinity, for any fixed e > there exists a > such that 

jx| > a =^ \a{x) — ao| < ^- (3-2) 

Cutting by additional Neumann boundary parallel to the y-axis at x = ±a, we get new 
operator Hj^^ defined using quadratic form. We can decompose this operator Hj^^ = 

Ha^c , where the "tail" part corresponds to the two halfstrips {\x\ > a) and the rest to the central 
part with the Neumann condition on the vertical boundary. Using Neumann bracketing, cf. [25\ 

Section XIIL15], we get H^^^ < in the sense of quadratic forms. 

We denote amin(a) := inf a{x). Since aess{H^^. t) = [-E'i("mm), oo) and H^^^. ^ < H^J 

\x\>a 

in the sense of quadratic forms, we get the following estimate of the bottom of the essential 
spectrum of the "tail" part 

^i(amm) < inf (Jess (^S^). (3.3) 

Since the spectrum of is purely discrete, cf. (TJ Chapter 7], the minimax principle gives 

the inequality 

inf aess(^?S^) < inf aess(i^a). (3.4) 

Since the assertion (j3.2p yields oq — e < amin and Ei is an increasing function of q, we have 

E'i(ao - e) < -E'i(amm)- (3.5) 

Giving together (jS.Sp . (j3.4p . and (|3.5p we get the relation Ei{ao — e) < inf (7css(-f^a)- The claim 
then follows from the fact that Ei = Ei{a) is a continuous function. ■ 

Lemma 4. // lim a{x) — oq = then [£^i(ao),oo) C acssiHa)- 

\x\~>oo 
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Proof. Let A € [i?i(ao), oo). We shall construct a sequence {ipn}'^=i satisfying the assump- 
tions of Lemma El We define the following family of functions 



where xi is the lowest transversal function and ipn{x) ■= n ^/'^ip{x/n — n) with satisfying 

1) G c^m, 

2) Vx G M, < ip{x) < 1, 

3) VxG(-l/4,l/4),v9 = l, 

4) Vx G M\ [-1/4,1/4], ip = 0, 

5) lly'llLZCM) = 1- 

Note that supp C (n^ — n, Ji^ + n). It is clear that tpn belongs to the form domain of Ha. 
The assumption (i) oi Lemma [2] is satisfied due to the normalization of xi and f. 

The point (ii) of Lemma [2] requires that {(j),'4jn) — > as n ^ oo for all (j) G C^(0), a dense 
subset of L?'{Q). However, it follows at once because (j) find ijjn have disjoint supports for n large 
enough. 

Hence, it remains to check that \\{Ha — A)V'nllD(h)* ^ as n ^ oo. An explicit calculation 
using integration by parts, boundary conditions (jl.2p and the relations 



n \m\L^{ 



(3.6) 



yields 



(a;,y)xi(y;ao)v3„(x)exp [i^JX- Ei{aQ)x) dxdy 



\{(t>,{Ha-X)ll,n)\ = 

+ 2i I (j){x,y)xi{y;cxo)ipnix)y/X- £'i(Qo)exp (za/ A - Ei{ao)x) dxdy 



< ^ll0llL2(n)llxi||L2({o,d))ll¥'llL2(K) + 2n ^\\(t)\/ X-Ei{ao)\\L2(^f^)\\xi\\Lmo,d))\\^\\L^R) 

for all (j) G D{ha)- The claim then follows from the fact that both terms at the r.h.s. go to zero 
as n ^ oo. ■ 



3.3 The existence of bound states 

Now we show that if a — oq vanishes at infinity, some behavior of the function a may produce a 
non-trivial spectrum below the energy £^i(ao). Note that this together with the assumption of 
Theorem [2] implies that the spectrum below £^i(ao) consists of isolated eigenvalues of finite mul- 
tiplicity, i.e. (TA\sc{Ha) 7^ 0. Sufficient condition that pushes the spectrum of the Hamiltonian 
below Ei{a{)) is introduced in the following theorem. 

Theorem 3. Suppose 

1) a{x) - ao G Li(M), 

2) /jg (q(x) - Qo) dx < 0. 

Then inia{Ha) < £'i(qo). 
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Proof. The proof is based on the variational strategy of finding a trial function ^lJ from the 
form domain of such that 

:= h^[i;] - Ei{ao)\ml2(^n) < 0- (3-7) 
Inspired by [18], we define a sequence 

i/jnix, y) := Lpn{x)xi{y] tto), Vn{x) := tT^I'^ ip{x j n) , 

where the function was defined in the proof of Lemma [H Using the relations (j3.6p and the 
integration by parts we get 

QalV'n] = f^~^llvi|lL2m-| + (|xi(0;ao)P + |xi(rf;"o)P) / ("(a^) - ao)v?„(x)dx. 

Since the integrand is dominated by the L^-norm of a — ao we have the limit 




by dominated convergence theorem. This expression is negative according to the assumptions. 
Now, it is enough to take n sufficiently large to satisfy inequality (j3.7p . ■ 

4 A 'rectangular well' example 

To illustrate the above results and to understand the behavior of the spectrum of the Hamiltonian 
in more detail, we shall now investigate an example. Inspired by [Hj we choose the function a to 
be a steplike function which makes it possible to solve the corresponding Schrodinger equation 
numerically by employing the mode-matching method. The simplest non-trivial case concerns 
a 'rectangular well' of a width 2a, 

I ai if < «, 
gl[x) = < 

I ao if |x| > a 

with a > and < ai < oq- In view of Theorem [3] this waveguide system has bound states. 
In particular, one expects that in the case when a\ is close to zero and oq is large the spectral 
properties will be similar to those of the situation studied in [9] . 

Since the system is symmetric with respect to the y-axis, we can restrict ourselves to the part 
of $7 in the first quadrant and we may consider separately the symmetric and antisymmetric 
solutions, i.e. to analyze the halfstrip with the Neumann or Dirichlet boundary condition at the 
segment (0, d) on y-axis, respectively. 

4.1 Preliminaries 

Theorem [2] enables us to localize the essential spectrum in the present situation, i.e., cJess 
[Ei{aQ),<x>). Moreover, according to the minimax principle we know that isolated eigenvalues 
of Ha are squeezed between those of H^i^c and Ha^\ the Hamiltonians in the central part with 
Neumann or Dirichlet condition on the vertical boundary, respectively. The Neumann estimate 
tells us that mia{Ha.) > -Ei(ai). One finds that the n-th eigenvalue En of Ha is estimated by 
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a. = 0.1 a, = 1 

1 1 




Figure 2. Ground state eigenfunctions for a/d = 0.3, ao = 20. 



4.2 Mode-matching method 

Let us pass to the mode-matching method. A natural Ansatz for the solution of an energy 
A e {Ei{ai),Ei{ao)) is 

OO ( COshjlnX) \ 

i^s/a{x, y) = 2_^ < \ , . > Xn(y; ai) for O < x < a, 



n=l 

OO 



sinh(Z„x) 
sinh((nii) 



^sla{^^ y) = X] exp(-/En(x - a))xn{y; ao) for x > 



n=l 



where the subscripts and superscripts s, a distinguish the symmetric and antisymmetric case, 
respectively. The longitudinal momenta are defined by 



In ■■= \/En{ai) - A, kn := y/Eniao) - A. 

As an element of the domain the function V should be continuous together with its 

normal derivative at the segment dividing the two regions, x = a. Using the orthonormality 
of {Xn} we get from the requirement of continuity 

OO r-d 

y^On / Xn{y,Oii)Xm{y;OiQ)Ay = hm- (4.1) 
n=l 

In the same way, the normal-derivative continuity at x = a yields 
°° f tanh 1 

^aJn^ coth 1*^^"") Xn{y;ai)xm{y;ao)dy + bmkm = 0. (4.2) 
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Substituting (j4.1|) to (j4.2p we can write the equation as 

Ca = 0, (4.3) 

where 



tanh 



/ Xn{y;ai)xm{y;ao)dy. 



d 







In this way we have transformed a partial-differential-equation problem to a solution of an 
infinite system of linear equations. The latter will be solved numerically by using in (|4.3|) an N 
by N subblock of C with large N . 



4.3 Numerical results 

The dependence of the ground state wavefunction with respect to ai for the central part of the 
halfwidth a/d = 0.3 is illustrated in Fig. [2l The slope of the wavefunctions at the boundary be- 
tween different Oj becomes smoother with increasing ai, at ai approaching oq the wavefunction 
is spread further and further in the x-direction, and at ai = qq bound state disappears. 

Fig. [3] shows the bound-state energies as functions of the 'window' halfwidth a/d in the 
case ao = 10^, ai = 10"^ (on the left) and ao = 20, ai = 5 (on the right). We see that 
the energies (in both cases) decrease monotonously with the increasing 'window' width and in 
the limit approach the uniform value Ei{ai). Their number increases as a function of a/d. 
Comparing this figure with Fig. 4 in [22] or Fig. 5 in |23], where the authors study the planar 
waveguide comprised of two straight regions connected by a region of constant curvature (with 
Dirichlet or mixed boundary conditions, respectively), we may see the similar behavior of the 
eigenvalues below the fundamental propagation constant. This comparison shows a community 
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0.05 



0.1 



0.15 
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0.25 



0.3 



a/d 



Figure 4. The dependence of first eigenvalue (in tlic units of (Tr/d)^) on a/d in the cases: (a) ag = 



50, ai = 3, (b) ao = 70, ai = 2, (c) 
Neumann case (e) adopted from [9]. 



100, ai = 1, (d) Qfo = 200, ai = 0.5, and in the Dirichlet- 



of the physical processes taking place in different structures where the longitudinal uniformity 
of the waveguide is broken by some obstacle - either by the different boundary conditions or the 
bend. For illustration, the first gap, i.e., the difference between first and second eigenvalue (or 
between first eigenvalue and the bottom of the essential spectrum if there is only one eigenvalue) , 
is plotted in the bottom part of Fig. [3l It is interesting that the dependence on a/d in the region 
where more than one eigenvalue exists is not monotonous, but there is a maximum. 

We can compare the dependence of the first eigenvalue on a/d with the system studied in [9] . 
The authors there are interested in the model of straight quantum waveguide with combined 
Dirichlet and Neumann boundary conditions. Fig. [4] shows that the larger oq and the smaller ai 
is (curves a, 6, c, and d), the closer to the Dirichlet-Neumann case (curve e) is the behavior of 
the first eigenvalue as a function of a/d. 
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